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We explore the robustness of universal dynamical scaling behavior in a quantum system near 
criticality with respect to initialization in a large class of states with finite energy. By focusing on 
a homogeneous XY quantum spin chain in a transverse field, we characterize the non-equilibrium 
response under adiabatic and sudden quench processes originating from a pure as well as a mixed 
excited initial state, and involving either a regular quantum critical or a multicritical point. We 
find that the critical exponents of the ground-state quantum phase transition can be encoded in the 
dynamical scaling exponents despite the finite energy of the initial state. In particular, we identify 
conditions on the initial distribution of quasi-particle excitation which ensure Kibble-Zurek scaling to 
persist. The emergence of effective thermal equilibrium behavior following a sudden quench towards 
criticality is also investigated, with focus on the long-time expectation value of the quasi-particle 
number operator. Despite the integrability of the XY model, this observable is found to behave 
thermally in quenches to a regular quantum critical point, provided that the system is initially 
prepared at sufficiently high temperature. However, a similar thermalization behavior fails to occur 
in quenches towards a multi-critical point. We argue that the observed lack of thermalization 
originates in this case in the asymmetry of the impulse region that is also responsible for anomalous 
multicritical dynamical scaling. 
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I. INTRODUCTION 

Characterizing the non-equilibrium dynamics of quan- 
tum many-body systems is of central significance to 
both condensed-matter physics and quantum statisti- 
cal mechanics. A quantitative understanding of non- 
equilibrium quantum phase transitions (QPTs) is, in par- 
ticular, a fundamental prerequisite for uncovering and 
controlling quantum phases of matter 1,2 , as well as for 
assessing the complexity of quantum annealing or adi- 
abatic algorithms 3 ' 4 . Unlike standard phase transitions 
which are driven by a change in temperature, QPTs are 
driven entirely by quantum fluctuations at zero temper- 
ature. They nevertheless share with their classical coun- 
terpart the generic feature of universality: in equilibrium, 
the critical properties of a system sufficiently close to a 
quantum critical point (QCP) depend only on a few rele- 
vant characteristics such as its symmetry and dimension- 
ality, thus defining the universality class to which the 
corresponding (continuous) QPT belongs. The univer- 
sality class is distinguished by a small set of critical ex- 
ponents - for instance, v and z, describing the power-law 
divergence of the characteristic length scale and the van- 
ishing of the characteristic energy scale, respectively 1 . 
In a non-equilibrium scenario, the system can be driven 
across a QCP dynamically, that is, through an explicit 
time-dependence of one or more control parameters in 
the underlying many-body Hamiltonian. This naturally 
prompts a number of questions: to what extent can uni- 
versal quantum scaling laws persist out of equilibrium 
and be solely specified in terms of the equilibrium phase 
diagram? Conversely, how does quantum criticality influ- 



ence the ability of a system to relax back to equilibrium? 

Historically, the first theoretical studies in these direc- 
tions date back to the pioneering work by Barouch and 
coworkers 5-7 which led, in particular, to the discovery 
of "non-ergodic" behavior in the zero-temperature long- 
time magnetization of a driven XY spin chain 6 . In re- 
cent years, the demand for quantitatively addressing the 
above broad questions has heightened under the impe- 
tus of experimental advances in manipulating ultracold 
atomic gases, which are enabling the unitary dynamics of 
many-body quantum systems to be explored with an un- 
precedented level of coherent control and isolation from 
the environment 8 ' 9 . As a result, non-equilibrium quan- 
tum critical physics is being actively investigated both 
from a theoretical and experimental standpoints. 

In this framework, an important step forward is pro- 
vided by the prediction of universal behavior in adia- 
batic dynamics near a QCP based on the Kibble-Zurek 
scaling (KZS) argument 10 (see also Ref. 11 for related 
independent work and Ref. 12 for a recent review). 
Originally introduced in the context of classical (finite- 
temperature) phase transitions in cosmology 13 , the KZ 
argument rests on the basic intuition that, irrespective of 
how slowly a system is driven across a continuous phase 
transition, adiabaticity is necessarily lost in the thermo- 
dynamic limit due to the vanishing energy gap at the 
critical point. Qualitatively, this determines typical time 
and length scales, t and £ respectively, that character- 
ize the adiabatic-to-diabatic crossover and, since "order" 
cannot be established on distances larger than £, results 
in the formation of a domain structure and the generation 
of a finite density of "topological defects" in the system. 
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Quantitatively, let the time-dependence in the quantum- 
mechanical Hamiltonian H(t) be introduced through a 
control parameter X(t), with A c = X(t c ) corresponding to 
the crossing of an isolated QCP at time t c , which can be 
taken to be t c = without loss of generality. If the sys- 
tem is initially (t = to) in the ground state, its ability to 
adiabatically adjust to H (t) is determined by the condi- 
tion that the typical time scale r(t) = \(X(t) — A c )/A(t)| 
associated with the applied control be sufficiently long 
relative to the slowest response time r r ~ A -1 , which is 
set by the smallest energy gap A. Since, as the QCP is 
approached, the latter vanishes as A ~ |A(f) — X c \ vz , 
adiabaticity is broken throughout an "impulse region" 
[t c — i, t c + t\ symmetrically located around the QCP, 
where the "freeze-out" time t is determined by the condi- 
tion t(£) = T r (t). In the simplest case of a linear sweep 
across the QCP (Fig. 1, top), A(i) — A c = tjr for a fixed 
rate t -1 > 0, this yields t ~ t vz ik vz + 1 ) and a typical 
gap A ~ £~ x . Correspondingly, the typical correlation 
length £ ~ ~ A -2 also scales with the quench time 
r. Since £ is the universal length scale near criticality, it 
determines the scaling of the final (t = tf) density of de- 
fects and, more generally, the total density of excitations, 
n ex (i/), created in the system. If d denotes the spatial 
dimension, the KZS result then follows: 



<(*/) 



_-dv/(vx+l) 
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The validity as well as the limitations of the above KZS 
have been carefully scrutinized in a number of settings. 
By now, the original KZS has been confirmed for a vari- 
ety of models involving a regular isolated QCP 14-18 , and 
extensions have been introduced for more general adia- 
batic dynamics, including repeated 19 , non-linear 20 , and 
optimal 21 quench processes. In parallel, departures from 
the KZS predictions have emerged for more complex adi- 
abatic scenarios, involving for instance quenches across 
either an isolated multicritical point (MCP) 20, 22-25 or 
non-isolated QCPs (that is, critical regions ) 26-29 , as well 
as QPTs in infinitely-coordinated 30 , disordered 31 , and/or 
spatially inhomogeneous systems 12 ' 32 . A main message 
that has emerged from the above studies is that, unlike in 
the standard KZS of Eq. (1) where the non-equilibrium 
critical exponent is completely specified in terms of static 
exponents, additional details about the time-dependent 
excitation may play an essential role in general. As a re- 
sult, genuinely non-static, path- dependent exponents may 
be required for dynamical scaling predictions. This fea- 
ture is vividly exemplified, for instance, in multicritical 
quantum quenches, whereby the asymmetry of the KZ 
impulse region relative to the static QCP (Fig. 1, bot- 
tom) causes a path-dependent minimum gap other than 
the critical gap to be relevant and an effective dynamical 
exponent Z2 ^ z to emerge 23 . 

In addition to characterizing the response to an adia- 
batic probe, the opposite limit of a sudden change of the 
tuning parameter X(t) near a QCP has also attracted a 
growing attention recently, in connection with the study 
of both dynamical quantum-critical properties 33-35 , as 



Adiabatic 




Adiabatic 



t r =0 



Adiabatic Impulse 




Adiabatic 
► 



~U t r =0 



t'2 



FIG. 1: (Color online) Qualitative sketch of the adiabatic- 
impulse-adiabatic sequence of regimes relevant to dynamical 
scaling arguments. Top: Symmetric impulse region, as as- 
sumed by the standard KZS scenario. Bottom: Asymmetric 
impulse region, as resulting from the existence of quasicritical 
path-dependent energy states, see Ref. 23. 



well as thermalization dynamics in closed quantum sys- 
tems and its interplay with quantum integrability . 
In particular, for sudden quenches with a sufficiently 
small amplitude, the existence of universal scaling behav- 
ior has been established for various physical observables 
and qualitatively related to the above KZ argument 35,44 : 
by associating to a quench with amplitude SX the charac- 
teristic length scale £ ~ |(SA| _l/ , and by interpreting £ as 
the correlation length in the final state, one immediately 
infers that n ex (t f ) ~ <r d - |<5A| dl/ . 

With a few exceptions where quenches at finite tem- 
perature and the associated thermal corrections have also 
been examined 35 ' 36,44 ' 45 , the large majority of the ex- 
isting investigations have focused on quench dynamics 
originating from the ground state of the initial Hamilto- 
nian H(to). Our goal in this work is to present a ded- 
icated analysis of finite- energy quantum quenches, with 
a twofold motivation in mind. Conceptually, elucidat- 
ing to what extent and how universal scaling properties 
may depend upon the details of the system's initializa- 
tion is needed to gain a more complete picture of non- 
equilibrium quantum-critical physics. While one might, 
for instance, naively expect that a sizable overlap with 
the initial ground state would be essential in determining 
the applicability of ground-state scaling results, a main 
highlight of our analysis is that the support of the initial 
state on those excitations relevant to the path-dependent 
excitation process is key in a dynamical scenario, in a 
sense to be made precise later. Furthermore, from a prac- 
tical standpoint, perfect initialization of a many-body 
Hamiltonian in its exact ground state is both NP-hard 
in general 46,47 and experimentally unfeasible due to lim- 
ited control. In this sense, our investigation both ex- 
tends previous studies on finite-temperature signatures 
of static QPTs 48 , and may be directly relevant to experi- 
ments using ultracold atoms 8 as well as nuclear magnetic 
resonance (NMR) quantum simulators 49 ' 50 . 

While our analysis focuses on the simplest yet paradig- 
matic case of an exactly solvable XY quantum spin chain 
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(Sec. II), we address non-equilibrium dynamics originat- 
ing from a large class of (bulk) initial states for a vari- 
ety of different quench schemes involving either a regular 
QCP or a MCP. Both pure and mixed initial states car- 
rying finite excitation energy above the ground state are 
examined, under the main assumption that, subsequent 
to initialization, the system can be treated as (nearly) 
isolated, hence evolving under a time-dependent Hamil- 
tonian. In particular, dynamical scaling in adiabatic and 
sudden quenches starting from an excited energy eigen- 
state are analyzed in Sec. III. A and III.B respectively, 
with emphasis on making contact with previously intro- 
duced adiabatic renormalization approaches 15,28 and on 
clarifying formal connections between scaling behavior in 
sudden and adiabatic dynamics. The case of a generic ex- 
cited pure state prepared by a sudden parameter quench 
is also considered in Sec. III.C, and criteria are identi- 
fied for KZS to be obeyed. Sec. IV is devoted to quench 
dynamics resulting from an initial thermal mixture, with 
the main goals of characterizing the robustness of dy- 
namical scaling behavior in realistic finite-temperature 
conditions, and of further exploring the conditions lead- 
ing to effective thermalization of certain physical observ- 
ables following a sudden quench toward criticality. In the 
process, we continue and extend the analysis undertaken 
in Deng et a/. 23 , by presenting finite-temperature gen- 
eralizations of the scaling predictions obtained for adia- 
batic (both linear and non-linear) multicritical quantum 
quenches, as well as evidence of how the peculiar nature 
of a MCP may also result in anomalous thermalization 
behavior. Section V conludes with a discussion of the 
main findings and their implications, along with further 
open problems. 



II. MODEL HAMILTONIAN 

A. Energy spectrum and equilibrium phase 
diagram 

We consider the homogeneous one-dimensional spin- 
1/2 XY model described by the Hamiltonian: 

N 

tf = -E(-r^ +1 + V^ +1 "^)' (2) 

where periodic boundary conditions are assumed, that 
is, a J a = o- J a +N , and N is taken to be even. Here, 
7, h G [—oo, oo], parameterize the degree of anisotropy in 
the XY plane, and the uniform magnetic field strength, 
respectively, in suitable units. The diagonalization of 
the Hamiltonian (2) is well-known 5,6 ' 51,52 , and we only 
recall the basic steps here. Upon introducing canonical 
fermionic operators {cj, cA via the Jordan- Wigner map- 



ping Cj = Ilf=i(~ cr z) cr +' H rewrites as a quadratic form 
N-i N 

h = -J2 ( c I c j+i + ^A+i + h - c -) + 2h J2 c h 

3 = 1 j=l 

- hN + V (4ci + 7 cjv ct i + h-c), (3) 

where the last term originates from the spin peri- 
odic boundary conditions and the parity operator V = 

Il3Li(-oi) = e in ^ c hi = +1(-1) depending on 
whether the eigenvalue of the total fermionic number op- 
erator is even (odd), respectively. Physically, V corre- 
sponds to a global Z2-symmetry which, for finite N, al- 
lows the even and odd subspaces to be exactly decoupled, 
H = H( + } + H^~\ and the diagonalization to be carried 
out separately in each sector. 

In finite systems, the ground state as well as excited 
energy eigenstates with an even number of fermions be- 
long to the V = +1 sector. By using a Fourier trans- 
formation to momentum space, c\, = x e~^ J cj, 
followed by a Bogoliubov rotation to fermionic quasipar- 
ticles {-fk,jl}, with 7fc = u k c k - iv k c'_ k , u k = u_ fe , 
V k = —V-kj ancl u k + v k = 1j the Hamiltonian in Eq. 
(3) rewrites as a sum of non-interacting terms: 

H(+) ^J2 Hk = E efe (^)^ + ^- fc - w 

k£K + keK + 

Here, the set K = A' + + A'_ of allowed momentum modes 
is determined by the anti-periodic boundary conditions 
on the fermions in the even sector, Cj+jv = — cj, which 

yields K± = { ± § ,±f , . . . , ±(tt - f ) }, and 

e k (h,j) =2v/(/i-cosfc) 2 + 7 2 sinfc 2 (5) 

is the quasi-particle energy of mode k. For each k, let 
H k = span{|O fc ), |l fe )}, where {|O fc ), |l fc ) = %\O k )} are or- 
thonormal states corresponding, respectively, to zero and 
one Bogoliubov quasiparticle with momentum k, that is, 

{Ok\ltlk\Ok) = 0, (lfc|7fc7fe|lfe) = 1, and similarly for 
—k. Thus, the four eigenstates of Hk provide a basis for 

B k = {|0 fe ,0_ fe )> life, 1-k), |0 fe , 1-fe), |l fe ,0_ fe >} (6) 



where the corresponding eigenenergies are given by — e^, 
efe, 0, 0, and a further separation into even (odd) sector 
for each k is possible due to the fact that \P k ,Hk\ = 0, 
with Vk = e irr ^ klk+1 - k ' 1 - k ^ = e l7Ti - c l Ck+c -k c -k) _ 

The ground state of H ^ corresponds to the BCS state 
with no Bogoliubov quasiparticles, 

l*o +) )= <S> |0fe0_ fc > = (g) (u fc + iv k cld k )\vac), 

k£K + keK + 

where |vac) is the fermionic vacuum. Many-body excited 
states in the even sector can be obtained by applying 
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pairs of Bogoliubov quasiparticle operators to |\&q + ^). In 
particular, excited eigenstates with support only on the 
even sector for each mode are obtained by exciting 
only pairs of quasiparticles with opposite momentum and 
have the form 

= ( ® |l*l-fe>)( (g) |0 fe 0_ fe >), (7) 

where labels the subset of excited modes. 

For finite N, the ground state and excited energy eigen- 
states with an odd number of fermions belong to the sec- 
tor V — — 1, which implies periodic boundary conditions 
on the fermions, c^+jv = Cj, and a different set K of 
allowed momentum modes, K = K + + K_ + {0, — it}, 

where K± = j ± ^, . . . , ±(tt - ^) j. Since one 

may show that ek=a = h — 2 and ek=—n = h + 2, occu- 
pying mode has always lower energy than occupying 
mode — 7r, thus the ground state of H is now 

lO = |l 0-*> (g) |0 fc 0_ fc ), 

and excitations may be generated by applying Bogoli- 
ubov quasiparticle operators in such a way that the con- 
straint on the total fermionic number is obeyed. Similar 
to modes in the even sector, k € K + , the subspace of 
each mode k £ K + yields four eigenstates of and, 
in principle, a basis formally identical to the one in Eq. 
(6) for the odd Hilbert-space sector. Although for finite 
N (thus necessarily in numerical simulations) V is al- 
ways a good quantum number under dynamics induced 
by H, the error in the computation of observables arising 
from identifying the two sets of modes K and K scales 
like 1/N. Thus, for sufficiently large N, a simplified de- 
scription in terms of a unique set of momentum modes is 
possible by using the basis 

B= (g) B k , (8) 

fee if + 

to characterize arbitrary states in the full Hilbert space 
H = &) keK+ {T-Lk ® H-k)- This becomes accurate in 
the thermodynamic limit N — > oo, where the many- 
body ground state becomes twofold degenerate and the 
Z2-symmetry spontaneously breaks, causing different V- 
sectors to mix. 

The equilibrium phase diagram of the model Hamil- 
tonian in Eq. (2) is determined by the behavior of the 
excitation gap of each mode, Afc(/i, 7) = ek{h,"f), with 
efe(/i,7) given in Eq. (5), and is depicted in Fig. 2. 
Throughout this work, we will mainly investigate scaling 
behavior in quenches involving either the regular QCP 
A (h c = 1,7 C = 1), which has equilibrium critical expo- 
nents v = z = 1 and belongs to the d = 2 Ising uni- 
versality class, or the MCP B (h c = 1,J C = 0), which 
has v = 1/2, z = 2 and belongs instead to the Lifshitz 
universality class 23 . In what follows, we shall refer to the 
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FIG. 2: (Color online) The phase diagram of the XY Hamil- 
tonian in Eq. (2). The regular QCP A (h c = l,j c = 1) and 
the MCP B (h c = 1, 7c = 0) are marked. The dashed-dotted 
(green) line separates the ferromagnetic (FM) and paramag- 
netic (PM) phases, whereas the dotted (blue) line represents 
the superfluid phase (SF). 

critical mode k c as the mode whose gap A&(/i, 7) vanishes 
in the thermodynamic limit. For both the QCPs A and 
B of interest, we thus have k c = in the large- N limit. 



B. Dynamical response indicators 

If the system described by Eq. (2) is driven across a 
QCP by making one (or both) of the control parameter(s) 
h, 7 explicitly time-dependent, excitations are induced as 
a result of the non-equilibrium dynamics. Since the gap 
vanishes at the QCP in the thermodynamic limit, this 
happens no matter how slow the Hamiltonian changes 
with time. In Refs. 17,28, the excess expectation value 
relative to the instantaneous ground state was shown to 
successfully characterize dynamical scaling behavior for a 
large class of physical observables in adiabatic quenches 
originating from the ground state. That is, for an exten- 
sive observable O, the following quantity quantifies the 
underlying adiabaticity loss: 

AO(t) ee (*(t)|0|*(t)>-(*(t)|0|*(t)), (9) 

where |\&(t)) and ^(i)) are the actual time-evolved state 
and the adiabatically evolved state resulting from |^(to)), 
respectively. For a generic quench process, where in prin- 
ciple both the time-dependence in H(t) and the initial 
state p(to) can be arbitrary, it is desirable to character- 
ize the response of the system in such a way that no 
excitation is generated by purely adiabatic dynamics 15 
and zero-energy quenches are included as a special case. 
This motivates extending the definition of Eq. (9) to 

AO(t) ee Tr[0(t)p(t)] - Tv[0(t)p(t)} , (10) 
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where now p(t) and p{t) are the actual time-evolved den- 
sity operator and the density operator resulting from adi- 
abatic evolution of p(to) respectively, and we also al- 
low, in general, for the observable O to be explicitly 
time-dependent. Let H(t)\^,(t)} = Ei(t)\^f i(t)) define 
snapshot eigenstates and eigenvalues of H(t) along a 
given control path. Then the adiabatically steered state 

P(t) = Ynj P<j(*o)l*i(*)><*j(*)|i with PiA*o) bein S the 
matrix elements of the initial state p(to) in the eigenstate 
basis \i£ i(to)) (ty j (ta)\ of the initial Hamiltonian H(to). 

With respect to the basis given in Eq. (8), a generic 
uncorrelated state in momentum space may be expressed 
in the form p(t) = &) keK+ Pk(t), where p k (t) is the four- 
dimensional density operator for mode k. Relative to a 
snapshot eigenbasis 

B k (t) = {\i/i(t))}, j = 0,... ,3, 

similar to the one given in Eq. (6), but constructed 
from time-dependent quasiparticle operators such that 

7fc(*)|0fc(*)> = 0, 7fe(*)|0fc(t)> = Mt)), Pk(t) may be ex- 
pressed as: 

i,j=0,3 

Suppose that the time-evolution operator for mode k 
is U k (t), that is, Pk (t) = U k {t) Pk {t )ul{t). Direct 
calculation shows that \0 k (t),l- k (t)) = c^ fc |vac), and 

|lfc(i), 0-k(t)) = cjjvac) for all t, which indicates 
that the snapshot eigenstates belonging to the Vk — 
— 1 eigenvalues are frozen in time, |0fc(t), l- k (t)) = 

|0 fc (*o),l-fc(to)), |lfc(t),0_ fc (t)> = |lfc(to),0_ fc (*o)). As 
long as Vk is conserved under H k (t), the even and odd 
sectors for each k are decoupled. Thus, upon letting 

u£(t)\U(t),l- k (t)) = a , k (t)\0 k (t o ),0. k (t )) 
+ a lik (t)\l k (t ),l- k (t )), 

we can evaluate the time-dependent excitation probabil- 
ity for mode k as follows: 

p k {t) = Tt[p k (tH(th k (t)} 

= Tir[p fc (t)(|l*(t) > l_ fc (t))(l fc (t),l_ fc (t)|] 

+ Tk[^ J fc(t)|l J fc(t),0-fc(t)}(lfc(*) > 0_ fc (t)|] 

= (poo,fc(*o) - Pn,k(ta))\ao,k(t)\ 2 + Pii,k(h) 

+ 2Re[p 0hk {t )a* ok (t)a hk (t)} + pas.fcfo), (H) 

where the relationships \ao, k (t)\ 2 + \ai yk (t)\ 2 = 1 
and pw.k = Pqi k have been exploited. Notice 
that from the above definition of ciQ_ k (t), it follows 
that |ao,/c(i)| 2 is the time-dependent probability that 
mode k is excited when it is in its ground state at 
t = to- Similarly, we may express the adiabatically 
evolved density operator p(t) = &) k <=K + Pk(t), with 

Pk(t) = T,ij=q,3Py,k(h)\i> k (t))(ijl(t)\. Thus, the time- 
dependent excitation probability of mode k relative to 



the adiabatic path is simply 

P k (t) = Tr[p k (t) 7 l(t) 7k (t)] 

= Pn,M + p33,fc(«o) = Pk(t ). (12) 

Upon combining Eqs. (11)-(12), the relative excitation 
probability of mode k is given by 

AP fc (i) = P k (t) - P k (t ) 

= (poo,fc(*o) - Pii,k(t ))\a ,k(t)\ 2 

+ 2Re[p Q1 , k (t )al k (t)a hk (t)]. (13) 

Physically, a non-zero contribution P k (to) may account 
for initial excitations due to either a coherent preparation 
into an excited state or to a finite temperature T. Clearly, 
if mode k is initially in its ground state, Pfc(to) = 0, we 
consistently recover the definitions in Ref. 28 for zero- 
energy quenches. Two relevant limiting cases of Eq. (13) 
will play a special role in what follows. First, if mode k 
is initially in a generic pure state of the form 

llM*o)> = J2 c Mih)), 

j=0,3 

then pm,k{h) = |co,fc| 2 , Poi,k(*o) = co,fccJ )fcI Pxi,fc(*o) = 
ci.fcj 2 , hence 

AP fe (f) = (|c , fc | 2 -| Cl , fc | 2 )|ao, fc (t)| 2 

+ 2Re[o, fc <$ iJfc a5 ifc (t)o llfc (t)]. (14) 

Second, if the initial state p(to) is a statistical mixture, 
then pio,k(h) = poi,k(<o) = 0, and we have 

AP k (t) = {pooM - Pn.k{to))\a .k{t)\ 2 ■ (15) 

The time-dependent excess expectation value AO(t) 
in Eq. (10) may be expressed directly in terms of the 
relative excitation probability for observables that obey 
[0(t),H(t)] = at all times. In this work, we shall pri- 
marily focus on the following choices: 

• 0{t) = Efc 6 «r+[7j[(t)7fc(*)+7lfc(*)7-fc(*)], leading 
to the relative total density of excitations: density: 

An cx (i) = f E ^[(Pk(t) - Pk(t)ht(th k (t)} 

k£K + 

which coincides with n ex (t) when the initial state is the 
ground state. 

• 0{t) = H(t), leading to the relative excitation energy 
density: 

AH ^ = | E n(Pk(t) - Pk(t))H k (t)} 

k£K + 

= 1 J2 e k (h(t)Mt))APk(t). (17) 

k£K + 
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While An ex (t) is especially attractive from a theory 
standpoint in view of its simplicity (possibly enabling 
analytical calculations), a potential advantage of AH(t) 
is that it may be more directly accessible in experiments. 

As a representative example of an observable not com- 
muting with the system's Hamiltonian, we shall addition- 
ally include results on the scaling behavior of: 

• O = XX = jj J2iLi <7 x <t x +1 j corresponding to the 
nearest-neighbor spin correlator per site along the in- 
direction 17 . In the Ising limit (7 = 1), the operator 
Af = (1 — XX) /2 is a natural measure for the "density 
of kinks" created by the quench, which directly relates 
to the number of quasi-particles excited at h = 12,14,40 . 
We have: 

AXX(t) = — A(-2cos&4c fc ) 



+ 



+ A(i7(t)sinfc(4cL fe -h.c.)) (18) 

In principle, the sums in Eqs. (16)— (18) should include 
all the modes in K +1 as indicated. However, to the pur- 
pose of analytically investigating dynamical scaling be- 
havior, it is useful to note that not all the allowed modes 
will necessarily change their state along the adiabatic 
quench path, effectively making no contribution to the 
relative expectation AO(t). In what follows, we shall re- 
fer to the subset of modes Kr C K + whose state changes 
in an adiabatic quench as the relevant modes. Let a 
power-law adiabatic quench process be parametrized as 
\{t) — A c = |i/r|"sign(t), a = 1 corresponding to the 
standard linear case also discussed in the Introduction. 
We may relate the number of relevant modes, Nr = | Kr \ , 
to the system size and the quench rate via Nr(N, t) oc 
-/V|fc max (r) — k c \, where fc max is the largest momentum in 
the relevant mode set. Since adiabaticity breaks at a time 
scale i ~ r «^/(«^+i) ) an d the typical gap, A ~ i' 1 , an 
accessible excited state contributes to the excitation if 
and only if its minimum gap along the path, A&, matches 
with this typical gap, A k ~ A. In general 23 , A& scales as 
Afe ~ (k — k c ) Z2 , where Z2 is a genuinely non-static expo- 
nent. Accordingly, the scaling of fc max can be determined 
by A - (fc max - k c ) Z2 , leading to 



z/[z 2 (avz + l)] 



(19) 



III. 



QUENCHES FROM A PURE EXCITED 
STATE 



A. Adiabatic quench dynamics from an excited 
energy eigenstate 

Adiabatic quenches from the ground state of the initial 
Hamiltonian H(to) have been extensively studied and are 
well understood in this model 14-17,23 . In order to explore 
the role of initialization, a first natural step is to investi- 
gate dynamical scaling behavior when the initial state is 
an excited eigenstate of H(to). Since, as remarked, the 



time-evolution of excited components along |0fe, 1-k) and 
|lfc, 0_fc) is trivial, we shall focus on excited energy eigen- 
states with support only on the even sector of each mode 
k, that is, on states of the form given in Eq. (7). Noting 
that there are only two possibilities for each mode, either 
c , fe = 1 or c ,a; = 0, Eq. (14) yields 

AP k (t) = (|c , fc | 2 - |c lifc | 2 )K fe (t)| 2 = ±|a , fe (t)| 2 , (20) 
and, correspondingly, 



A7i cx (t/) = ^ J2 [±l a o,fe(i/)P 



(21) 



Thus, the relative excitation density is the same, up to 
a sign, in two limiting cases: (i) the many-body ground 
state, corresponding to co.fc = 1 for all k and to an over- 
all positive sign in Eq. (21); and (ii) the state where all 
allowed pairs of quasiparticles are excited, correspond- 
ing to co,fc = for all k and to an overall negative sign 
in Eq. (21). Since KZS is known to hold for a linear 
quench process with initial condition (i), and a global 
sign difference would not change the scaling behavior, 
KZS is expected to persist for the maximally excited ini- 
tial eigenstate (ii) as well. This is to some extent surpris- 
ing both in view of the fact that such an initial state has 
zero overlap with the BCS state |\I/q + ' > ) , and because one 
would not a priori expect highly energetic eigenstates to 
be sensitive to the ground-state QPT. 

Interestingly, critical properties of excited eigenstates 
in the XY chain have recently attracted attention in 
the context of static QPTs 53 . Suppose that each ex- 
cited eigenstate is associated with an ordered binary 
strings of length \K\ = N, where (1) represents a 
mode in its ground (excited) state, respectively. Then 
a compact description of the eigenstate may be given in 
terms of the "discontinuities" of a (suitably regularized as 
N — > 00) characteristic function of the corresponding oc- 
cupied mode set, where no discontinuity is present when 
all modes are or 1, and a discontinuity is counted ev- 
ery time the occupation of a given mode changes along 
the string. Alba et a/. 53 have analytically proved, in par- 
ticular, that the block entanglement entropy, Sl, of an 
excited eigenstate of the critical XY chain may still obey 
conformal scaling as in the ground state provided that 
the number of discontinuities remains finite in the ther- 
modynamic limit, that is, Sl ~ logL, where L 1 
is the block size. Conversely, Sl exhibits non-critical 
scaling, Sl ~ L, when the number of discontinuities be- 
comes itself an extensive quantity as N — > 00. Thus, 
certain highly excited states (including the fully excited 
state considered above) can still display critical behav- 
ior, the number of discontinuities in the full set of mo- 
mentum modes being the essential factor in determin- 
ing the static scaling behavior. While, intuitively, non- 
analyticities in the characteristic mode occupation func- 
tion need not play a direct role for simple observables 
such as the excitation density, these results still prompt 
the following question: to what extent does a distinction 
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FIG. 3: (Color online) Scaling behavior of the final relative 
excitation density in a linear quench of the magnetic field h 
around the QCP A in the Ising chain, starting with different 
excited eigenstates of H(h c ). Panel (a): only k c — n/N is 
excited initially. The linear fit for N = 3200 yields -0.535 ± 
0.002. Panel (b): the five lowest-energy modes are initially 
excited. A linear fitting slope of — 0.549±0.003 is now reached 
at N = 12800. Panel (c): five modes (k = k c , 5tt/N, 9iy/N, 
13n/N, and 17n/N) are initially excited. The linear fit for 
N = 12800 yields -0.546 ± 0.002. Panel (d): the five lowest- 
energy modes are initially excited for N = 12800 as in (b), but 
as the system size is increased linearly, the number of excited 
modes is increased accordingly. In all cases, the relevant r- 
range r m i n < 20 < r < 250 < r max (see text and Fig. 4). 



between "critical" (leading to KZS) and "non-critical" ex- 
cited eigenstates exist for dynamical QPTs? 

A key difference with respect to the static situation 
is that only the relevant modes matter in a dynamical 
QPT, k e Kr = [k c , kjoax], with £; max given in Eq. (19). 
In Fig. 3, we present exact numerical results, obtained 
by direct numerical integration of the time-dependent 
Schrodinger equation, for the relative excitation den- 
sity in a linear adiabatic quench of the magnetic field 
h around the QCP A of the critical Ising chain (7 = 1). 
Different initial eigenstates are compared over a common 
range of r, which is chosen to be well within the appropri- 
ate r-range 28 for ground-state quenches (see next para- 
graph and Fig. 4 for further discussion of this point). 
In panel (a), the system is initialized in the first ex- 
cited state, where only the critical mode is initially ex- 
cited (thus only one discontinuity is present), whereas 
in panel (b), the five lowest-energy modes are initially 
excited (leading to one discontinuity as well). In case 
(a), while no scaling is visible for a system with size 
A" = 400, progressively better scaling behavior emerges 
as A" is increased, with the value at A" = 3200 approach- 
ing the asymptotic KZS value (and better agreement be- 
ing achievable by optimizing the r-range, see below). In 
constrast, for the data in panel (b), a system size as large 
as JV = 12800 is required for a scaling of comparable 
quality to be established. Since the only difference be- 



tween cases (a) and (b) is a different (fixed) number of 
initially excited modes in Nr, the fact that upon increas- 
ing A" (thereby increasing Nr accordingly) a better KZS 
is comparatively obtained in (a) suggests that the ratio 
between the number of initially excited (or non- excited) 
modes and Nr is crucial for dynamical scaling behavior 
- not (as intuitively expected) the discontinuity proper- 
ties which characterize the initial mode occupation per 
se. More explicitly, let Ne denote the number of modes 
in Kr that are excited at time t$, with Nr. — Ne corre- 
spondingly denoting the number of non-excited modes in 
Kr, and let 

Mr = mhx{N E , N R — Ne}- 

Motivated by the above observations and also recalling 
the symmetric role played by initially non-excited vs. ex- 
cited modes in determining the time-dependent relative 
probability of excitation [Eq. (20)], we conjecture that 
KZS emerges in the thermodynamic limit provided that 
the initial excited eigenstate satisfies 



Mr 
Nr 



£<1. 



(22) 



Clearly, the case of ground-state initialization corre- 
sponds to Ne = Mr = 0, and the fully excited state 
coincides with Ne = Nr, Mr = 0. For a generic initial 
excited eigenstate, Eq. (22) allows in principle Mr to be 
an extensive quantity in the thermodynamic limit. Two 
additional results are included in Fig. 3 to illustrate the 
above possibility. In panel (c), we still have five excited 
modes in Nr as in (b), but five discontinuities as opposed 
to just one. For the same system size (thus also the same 
e) , the scaling is not worse than in panel (b) , further sup- 
porting the conclusion that the number of discontinuities 
does not play a role towards the emergence of dynamical 
scaling. In panel (d), a fixed value of e, equal to the one 
in (b) at A" = 12800, is explored for different values of N, 
by also proportionally increasing Ne- As the data show, 
the resulting An ox is the same, indicating that Mr may 
indeed be allowed to be an extensive quantity as long as 
Eq. (22) is obeyed. 

It is important to address how the choice of a range 
of r-values affects the above scaling conclusions. Let 
Tmin < t < r max and f min < t < f max denote the valid 
range for ground-state initialization 28 , and for excited- 
state initialization, respectively. Since Tmin is determined 
from the requirement that an adiabatic regime exists 
away from criticality, whereas r max follows from ensur- 
ing that adiabaticity can be broken in a finite-size sys- 
tem, both Tmin an d T max are related to the scaling of the 
many-body gap between the ground state and first (avail- 
able) excited state. Thus, f m j n (T max ) could in general be 
substantially different from T m j n (r max ), respectively. In 
our case, however, the Hamiltonian in Eq. (2) can be ex- 
actly decoupled into two-level systems for each mode k. 
Therefore, the relevant gap is always Afc c = £fe c (/i, 7), ir- 
respective of the initial condition. For this reason, the 
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FIG. 4: (Color online) Scaling behavior of the final relative 
excitation density in a linear magnetic-field quench around 
the QCP A in the Ising chain, starting with an eigenstate of 
H(t c ) where the five, ten, and twenty lowest-energy modes 
are initially excited for N = 12800, N = 25600, N = 51200, 
respectively [same as in Fig. 3(d)]. The relevant r-range is 
now f m in ~ r min = 5 < r < f max = 20 < r max . A linear 
fitting slope of —0.5019 ± 0.002 is now reached for all these 
cases, in agreement with the KZS prediction. 



any finite-energy initial state might imply more restric- 
tive constraints as compared to the zero-energy case. In 
particular, according to Eq. (22), not all the excited 
eigenstates can lead to KZS, and the better Eq. (22) is 
satisfied, the closer KZS will be approached. This ex- 
plains why, for instance, the fitting slope —0.549 from 
panel (b) of Fig. 3 is not as close to the KZ value as the 
one obtained for a ground-state quench with the same 
range of r. In the setting of Fig. 3(b), the majority of 
the relevant modes stay in their ground state. In order 
to reduce the contribution to An ex from the five lowest- 
energy modes, we can decrease t such that Nr will be 
increased in Eq. (22). Numerical support for this strat- 
egy is shown in Fig. 4, where an optimal range of r is 
identified for the same initial states as in Fig. 3(d), and 
very good agreement with KZS is recovered. Thus, we 
can conjecture that if the majority of modes that enter 
Mr are low-energy modes, we can reduce their contribu- 
tion to An ex by decreasing the upper bound to r. That 
is, we choose r min < t < f max , where f max < r max . Con- 
versely, if the majority of modes that are enter Mr are 
from high-energy modes, then we reduce the contribu- 
tion from these modes by increasing the lower bound to 
t. That is, we let f min < r < r max , with f min > r min . 

Additional theoretical understanding of the criterion 
given in Eq. (22) may be sought by invoking the per- 
turbative Adiabatic Renormalization (AR) approach 54 , 
which was successfully applied to explain the scaling re- 
sults for adiabatic quenches starting from the ground 
state 15,28 . Can first-order AR still capture dynamical 



scaling for initial excited eigenstates? Let us focus on lin- 
ear quenches (a = 1), and let the time-dependent Hamil- 
tonian be parametrized as H(t) = H c + ]A(£) — AJizi = 
H c + (t/r)Hi, with H c quantum-critical in the thermo- 
dynamic limit, so that the relevant QCP is crossed at 
t c = 0. If the system is prepared in the £-th eigestate of 
H(t(j), with to — > t c as in the examples previously con- 
sidered, the time-evolved state from first-order AR may 
be expressed in the form 



\*W(t)) = e- iT *W\*t{t)) 



£<#(*)l**»(t)}, 



where Ygit) includes in general both the Berry phase 

and the dynamical phase, and {t) gives the time- 
dependent amplitude along the m-th snapshot eigenstate. 
Following a derivation similar to the one given in Refs. 
28,54, and letting A m (i) = E m {t) - E e (t), we find: 



(t) 



-ir m (t) 



, (* m (f ) | Hi | tf , (i' ) ) i jf dsAm ( S ) 

dt ; — ; ; — : t J *° . 

tQ E m (t>) - E t (f) 



Thus, to first order in the quench rate 1/r the 
adiabaticity loss can be quantified by AO(t) = 
(I'M (t) \Q(t) | W (*)) - (*t (t) 1 0(t)\* t (*)>• In particular, 
this yields 

An«(f) = ! £ |cW(*)[ 2 ((* m (*)| £ 7+(i)7*(t)[<M*)) 



keK, 



Since Hi is a one-body operator in our case, 
the only non-zero matrix elements {^> m (t)\Hi\^/g(t)) 
include many-body eigenstates \^ m (t)) which dif- 
fer from in the occupation of precisely 

one mode. Thus, (*m(*)| Efc 7jfe(*)7fc(*)l*m(*)> - 

(*/(*)! E*7fc(*bfc(*)l*/(f)> = ±!> which im P lies 



A «e X (t/) = ^E [±I C ^(*/)| 2 ' 



(23) 



Except for a possible sign difference for each m, the above 
expression is formally identical to the one holding for 
ground-state initialization {£ = 0), in analogy with the 
exact Eq. (21). Numerical calculations of the relative ex- 
citation density according to Eq. (23) [for instance with 
the same initial condition as in Fig. 3(a), data not shown] 
confirm that the condition for initial excited eigenstates 
to support KZS is the same in first-order AR as the one 
conjectured based on exact numerical results. 



B. Sudden quench dynamics from an excited 
energy eigenstate 

As mentioned in the Introduction, scaling results for 
sudden quenches of the control parameter A around its 
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FIG. 5: (Color online): Scaling behavior of the final excitation 
density in combined magnetic-field ground-state quenches 
across QCP A in the Ising chain. Top: Sudden quench 
h c i y hf (see text) followed by a linear quench back to h c , 
with the system finally kept at h c . Bottom: Linear quench 
from h c followed by a sudden quench hf i-> h c , with the sys- 
tem finally kept at h c . In both cases, TV = 400. 



critical value A c have been recently obtained by De 
Grandi et al. 35 under the assumptions that the system 
is in the ground state of the initial Hamiltonian and the 
quench has a small amplitude, leading to a final excita- 
tion density 



t (t/) ~ |A- A, 



5X 



du 



(24) 



with SX <C 1 in suitable units. Before addressing, in 
analogy to the case of adiabatic dynamics, the extent 
to which the expected scaling behavior may be robust 
against initialization in a finite-energy eigenstate, it is 
useful to explore more quantitatively the connection be- 
tween ground-state adiabatic vs. sudden quenches im- 
plied by Eq. (24). 

Suppose, specifically, that the amplitude of a sudden 
magnetic-field quench near the QCP A of the Ising chain 
is directly related to the rate r of a corresponding linear 
adiabatic sweep across the same QCP via hf = h c ± h, 
where h oc t/r and t is the KZ freeze-out time scale, i.e., 



t 



z+1 >. Equation (24) then yields 
n cx (t f ) ~ \h f 



I d\J 



dv/(vz+l) 



(25) 



In other words, the scaling behavior resulting from 
Eq. (24) is essentially equivalent to KZS. While this 
could be quantitatively demonstrated by direct cal- 
culation of n cx (t) in a sudden quench, it can also 
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FIG. 6: (Color online) Scaling behavior of the final relative 
excitation density in a sudden magnetic-field quench across 
QCP A in the Ising chain, starting with the first excited state 
of H(h c ). The linear fitting slope for N = 3200 is -0.5244 ± 
0.0004 for 20 < r < 250. Closer agreement with the KZS may 
be reached by optimizing over r as in Fig. 4. 



be nicely illustrated by examining combined sudden- 
adiabatic quenches, which have not been explicitly ad- 
dressed to our knowledge 55 , and will also be relevant in 
Sec. IIIC. Two possible "control loops" starting from 
h(t Q ) = h c are depicted in Fig. 5: we can either (i) 
suddenly change the magnetic field amplitude h c t- > hf, 
and then adiabatically change it back to h c (top panel); 
or we can (ii) slowly ramp up h c to hf, and then sud- 
denly quench hf i— > h c (bottom panel). As it is clear 
from the numerical data, the total excitation density cre- 
ated from the combined sudden-adiabatic quench shows 
KZS throughout the entire process in both cases, pro- 
vided that t is within the appropriate scaling range 
Tmin < t < T max . Notice that the quench process de- 
picted in Fig. 5 is similar to the repeated linear quench 
across QCP A studied in Ref. 19, in the sense that the 
initial and final value of the control parameter coincide. 
While KZS was found to hold in such a repeated lin- 
ear quench, the difference is now that half of the linear 
adiabatic sweep is replaced by a sudden quench. Since, 
however, the interval [h c —h, h c +h] corresponds to the im- 
pulse region in the KZ scenario for a pure linear quench, 
the scaling results of the combined quenches under con- 
sideration may be understood as a consequence of the 
fact that the sudden quench component can only further 
enforce the impulse mechanism by which excitation is 
generated in the KZS argument. Interestingly, as long as 
the scaling exists, we can also observe that (i) and (ii) 
lead to almost the same final excitation density, even if 
the intermediate values of the excitation density after the 
sole sudden [in (i)] or linear [in (ii)] quench are different. 
In summary, the existence of KZS in ground-state sud- 
den and combined sudden-adiabatic quenches with small 
amplitude is essentially a reflection of the fact that the 
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system goes through an impulse region around the QCP 
no matter how slow or fast the quench is effected. 

While sudden quenches of arbitrary amplitude will be 
further considered in the next Section, we now return to 
the question of whether dynamical scaling also holds in 
small-amplitude sudden quenches when the system is ini- 
tially prepared in an excited eigenstate of H(to) = H c . 
Exact numerical results are presented in Fig. 6, where in 
order to ease the comparison with a linear quench, we 
have again explicitly related the sudden-quench ampli- 
tude to r as h f - h c oc t" 1 / 2 . The data for N = 3200 
indicate that the scaling exponent is slightly closer to the 
KZS prediction than the one obtained in a pure linear 
quench with the same initial condition and r-range [cf. 
Fig. 3(a)]. Since a sudden quench effectively strength- 
ens the impulse mechanism in the KZS argument, the 
number of relevant modes Nr is larger than the one in- 
volved in an adiabatic linear quench. Thus, for the same 
initial condition (the same Mr), the ratio e in Eq. (22) 
is smaller in a sudden quench than in a linear quench, 
comparatively leading to a scaling exponent closer to 
KZS. Therefore, our conclusions for excited-state sud- 
den quenches are consistent with the ones reached for 
excited-state adiabatic quenches, and reaffirm how small- 
amplitude sudden quench dynamics and adiabatic dy- 
namics near a QCP are essentially equivalent over a wide 
range of initializations. 



C. Adiabatic dynamics following a sudden quench 
from the ground state 

In addition to eigenstates of the initial Hamiltonian, 
another physically relevant class of initial preparations 
is provided by pure states that are reachable from the 
many-body ground state via a sudden parameter quench 
of arbitrary amplitude. For concreteness, let us focus 
on adiabatic dynamics following a sudden quench of the 
magnetic field h to its critical value h c in the Ising chain. 
Thus, the initial state for the adiabatic quench is a super- 
position of different eigenstates of the Hamiltonian H(h c ) 
after the (instantaneous) sudden quench. Since for each 
mode k the parity quantum number V k is conserved, and 
the ground state of H k lies in the even sector Vk = 1, the 
expansion coefficients C2. k = c%,k = 0, whereas co >k and 
Ci t k are obtained from expanding the ground state before 

the sudden quench in the eigenstate basis {l^jU^o ))} °f 
the quenched Hamiltonian H(h c ). 

We can picture the resulting dynamics in terms of a 
combined sudden-adiabatic quench process (see Fig. 7, 
inset), except that unlike in Sec. IIIB we only focus 
on the scaling behavior of the relative excitation density 
An ex (t) created after the sudden quench. Exact numeri- 
cal results are plotted in the main panel of Fig. 7, showing 
that for a large range of sudden-quench initializations, 
the final excitation density still obeys the same KZS, 
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FIG. 7: (Color online) Scaling of the final relative excitation 
density in an adiabatic magnetic-field quench across QCP A 
in the Ising chain, starting from an excited state prepared 
by suddenly quenching ho M- h c for different initial values 
of ho- The combined control path is illustrated in the in- 
set. The linear fitting slope for ho = -1,0.2,0.75,0.85,0.95 
is -0.50283 ± 5.0 x 10" 5 , -0.50697 ± 6.0 x 10" 5 , -0.5237 ± 
1.0 x 10~ 4 , -0.52800 ±5.0 x 10~ 5 , and -0.5037 ±8.0 x 10~ 4 
respectively. In all cases, the system size N = 400. 

as in adiabatic dynamics starting from the ground state. 
The above scaling result can be derived analytically in 
two limiting cases, starting from Eq. (14). Upon inte- 
grating over all the relevant modes, we find 

An cx {t) = -[ Tp k (t)dk = f maX {(2|c , fc | 2 - l)|a , fc (i)| 2 
Wo Jo 1 

± 2Re[c , k cl k al k (t) ai . k (t)}\ — . (26) 

There are two contributions in AP k (t). If the initial state 
of mode k is close to either a non-excited or to a fully 
excited state ( |co,fe | 2 ~ 1 or co.fc| 2 ~ for all k € Kr, 
respectively), the first term is the dominant one. In this 
case, KZS clearly holds. In the opposite limit where each 
mode k £ Kr is initially half-excited ( | Co , fc | 2 ~ 1/2), 
the second term is the dominant one. Since, for a sudden 
quench to h c , the latter is the center of the impulse region 
(recall Fig. 1, top) and at most half of the impulse region 
can be crossed, all the relevant modes can at most be 
close to half-excitation, making this second limiting case 
directly relevant to the sudden-quench state preparation 
for suitable ho. Assuming that \co, k \ 2 ~ 1/2 and ignoring 
relative phases thus yields 

AP fe (i) ~ \ao, k {t) ai ,k{t)\ ~ |a ,fcW|^/l-|ao,fc(t)| 2 - 

By invoking the Landau-Zener formula 16 , the asymptotic 
(t / — > oo) excitation probability for modes near k c scales 
like e ~ 27rk T when t -t -oo. Starting from QCP A (the 
center of the impulse region) will not, however, affect the 
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exponential behavior 56 . Therefore, \ao,k(t)\ 2 ' 
as long as tf is deep in the adiabatic region, and 1 



|oo,fc(*)|' 



k t. Integrating over the relevant modes 



then gives the anticipated KZS result: 



J dk \aa,k{t)\sj (1 - |ao,fc(*)| 



-1/2 

dk k 



-1/2 



-1/2 



.-1/2 



[Eq. (19)] in 



where we used the fact that fc max 
the upper integration limit. 

While the above argument suffices to explain the emer- 
gence of KZS starting from special sudden-quench initial- 
izations, for generic quenches the dominant term in Eq. 
(26) need not be the same for different modes. In order 
to gain further insight, it is necessary to inspect the dis- 
tribution of the excitation probability for each relevant 
mode after a sudden quench from a generic value ho i-> 
h c . Numerical results for the low- lying modes are pre- 
sented in Fig. 8 for a wide range of initial magnetic-field 
strength h$. For each mode k, we can identify two bound- 
ary values, h™'™ and h™*, such that when < h < 



i c , mode k is close to its ground state after the sudden 



mm 
0.A- 

2 



or ho > h 
w 1/2). Since 



max 
0,k » 



quench ( | Co , fc | 2 ~ 1), whereas if ho <C h 
mode k is close to half-excitation (|co 
^o"" an< ^ ^cTfc* are approximately symmetric with respect 

let us for simplicity take 
um Qualitatively, h Q u k 



to the critical value h c = 1. 

Kk = h o^> with h tk~ 2h c-Kk- 

can be determined by the condition A(/iq, k) w A(h c , fc), 



which yields approximately co.fc| ~ 1. If, conversely, 
A(h ,k) > A(h c ,k) (\h - h c \ > \h r Q n k - h c \), we can 
consider | Co , fc | 2 ~ 1/2. Altogether, the results in Figs. 7- 
8 indicate that the limiting analytical condition of re- 
quiring the same dominant term in Eq. (26) for all the 
relevant modes is too strong for An cx (t) to show KZS. 
For instance, when ho = 0.95, not all the relevant modes 
are staying in their ground state (fc c is not), yet KZS 
holds. In general, however, the variation of | co,fc | 2 with 
k does affect the scaling result. For instance, when ho 
is around 0.75, agreement with the KZS prediction for 
the same system size is relatively poor, motivating one 
to roughly identify the range 0.6 < ho < 0.9 with a cross- 
over region. Based on these observations, we conjecture 
that a necessary (and sufficient) condition for the rela- 
tive excitation density An cx (t) to approach KZS in the 
thermodynamic limit is that the dominant term in Eq. 
(26) is the same for the majority of the relevant modes. 

An alternative physical interpretation of the above con- 
jecture may be obtained by observing that for a generic 
value of ho, there exist modes k e ,k g <E K+ such that if 
k c < k < k e , A(ho, k) > A(h c , fc), whilst if k g < fc < tt, 
A(/io,fc) ~ A(h c ,k), and we also assume fc e < k g for 
concreteness. Since, in an adiabatic sweep with speed r, 
the set of relevant modes Kr = [fc c , fc ma x] is determined 
according to Eq. (19), we can distinguish three different 
regimes depending on how fc max is positioned relative to 
the interval [fc e ,fc g ]: 

(i) fcmax < k e : in this case, all the relevant modes 
are half-excited, recovering one of the limiting situations 




FIG. 8: (Color online) Dependence of the excitation coeffi- 
cient 2|co,fc| 2 — 1 upon the initial magnetic-field strength ho 
in a state prepared by a sudden quench ho i-> h c in the Ising 
chain (7 = 1). The values co,fc are obtained by expanding the 
ground state of H(ho) in terms of the eigenbasis of H(h c — 1) 
at QCP A. The five lowest-energy modes are considered, for 
system size N = 400. 



(analytically) leading to KZS, as already discussed (e.g., 
h = -l in Fig. 7); 

(ii) fc e < 

fcmax < fcg : in this case, by a reasoning similar 
to the one leading to Eq. (22), KZS is predicted to emerge 
provided that (fc m ax — fce) = efcmax, e <C 1, in such a way 
that the majority of the relevant modes are half-excited 
(e.g., h = 0.2 in Fig. 7); 

(iii) kg < fcmax : in this case, KZS is predicted to emerge 
provided that k g = efc m ax, e <C 1, in such a way that the 
majority of the relevant modes stay in their ground state 
(e.g., ho = 0.95 in Fig. 7). 

Thus, for both ho = 0.75 and ho = 0.85, the initial 
state prepared by the sudden quench may be interpreted 
to lie in the cross-over region between cases (ii) and (iii), 
explaining why the resulting scaling deviates appreciably 
from the KZ prediction. 

Similarly to the excited-eigenstate initialization, 
sudden-quench initialization will also add more con- 
straints on the appropriate r range for KZS to hold. If 
the initial state is prepared via a sudden quench that 
guarantees one of the above conditions (i)-(iii) to be ful- 
filled for any r £ [T m j n , r max ], then the latter range is also 
appropriate for KZS to emerge under excited-eigenstate 
initialization. If not, the situation is more involved, and 
the range of t may need to be adjusted such that either 
(ii) or (iii) is enforced. If condition (ii) is more likely to be 
obeyed (e.g., if ho ~ 0.6), we can choose f m i n > r m j n in 
such a way that the number of modes between fc e and k g 
is decreased, and the majority of relevant modes is thus 
half-excited. If instead condition (iii) is more likely to be 
obeyed (e.g., if ho ~ 0.9), we can choose f max < T max in 
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such a way that the number of relevant modes staying 
in their ground state is increased. While the strategy for 
adjusting the r-range in a sudden-quench initialization 
is similar to the one advocated in excited-eigenstate ini- 
tialization, conditions (i)-(iii) are in fact easier to fulfill 
than Eq. (22). For instance, for N = 400, the worst scal- 
ing in Fig. 7 is still relatively close to KZS, whereas the 
latter is completely lost when initially only k c is excited 
in Fig. 3(a). This difference is due to the fact that the 
initial occupation of modes in the relevant set changes 
less abruptly in a sudden-quench initialization than in 
excited-eigenstate initialization. 

We conclude our discussion of quench processes orig- 
inating from a (pure) excited state by commenting on 
the fact that the analysis developed for An cx (t) can be 
extended to different observables without requiring ma- 
jor conceptual modifications. While an explicit example 
involving the spin correlator defined in Eq. (18) will be 
included in the next Section, the basic idea is to pro- 
ceed in analogy with ground-state quenches 28 , by taking 
into consideration the appropriate scaling exponent as 
determined by the physical dimension of the observable 
O. Consider, for instance, the relative excitation energy 
AH(t) defined in Eq. (17) which, as remarked, can be 
experimentally more accessible than the relative excita- 
tion density. In all the situations where KZS holds for 
the latter, An ex (i/) ~ r -<W("*+i) ^ T -i/2 ( m particu- 
lar, in the case of excited-state initialization via a sudden 
quench just discussed), we also find for our model that 

AH(t f ) ~ T -W+*W("*+i) ^ T -\ 

consistent with the corresponding ground-state scaling 
behavior explored in Refs. 23,28. 



IV. QUENCHES FROM A THERMAL STATE 

A. Adiabatic quench dynamics 

While we have only focused thus far on initialization 
mechanisms resulting in a pure excited state, another 
large class of initial states with a finite excitation energy 
may be obtained through dissipative means, in partic- 
ular because the system may find itself (or be placed) 
in contact with a thermal bath. After a time sufficient 
for equilibration to occur, the system would then relax 
to a canonical ensemble at temperature T. In equilib- 
rium, it is well known that the influence of a ground- 
state QCP can cross over to a finite range of tempera- 
tures, the so-called "quantum critical regime," which is 
often broader than naively expected 1,57 ' 58 . In a dynami- 
cal scenario, how robust is dynamical scaling (in partic- 
ular, KZS) to initialization at a finite-temperature? If 
scaling persists, how do the relevant non-equibrium ex- 
ponents depend upon the initial temperature? Motivated 
by these questions, scaling behavior in a system initially 
prepared in a thermal equilibrium state at criticality and 



then adiabatically quenched away from the QCP has been 
analyzed in Ref. 35. In particular, it is shown that for 
fermionic quasi-particles, the excess excitation due to a 
quench across a standard QCP obeys 

Anc X (t/)~^T-( d +^/0- +1 \ (27) 

provided that the initial temperature is high enough 
(T 3> ek(to), for all k 6 Kr). Our goal here is to both 
present quantitative evidence for the above scaling law 
and, most importantly, to extend the analysis to multi- 
critical QCPs. 

Let T denote the initial thermal equilibrium temper- 
ature, so that the initial density operator has the form 
p(to) = (8)fceA' + Pk{to), with pk(t ) given by: 

PooAto) = ±e+^/ T , p lhk (to) = ±e-^» T , 
P22,fe(io) = PssAto) = \, (28) 

in units where h — ks — 1 and with 

2 = 2 + e+ e k(hn)/T + e -e fc (fe, T )/T_ 

For clarity, we focus on linear adiabatic dynamics first. 
We shall study both the standard Ising QCP A under 
a magnetic-field quench of the form h(t) = 1 — t/r 
[h = h c = 1,7 = 1 in Eq. (28)], and the MCP B 
under a simultaneous quench of the magnetic field and 
the anisotropic parameter, h(t) = 1 — j(t) = 1 — t/r 
[h = h c = 1,7 = 7 C = 1 in Eq. (28)]. At T = 0, 
the scaling of the excitation density can be in both cases 
described by n ex (*/) ~ T -duz/[z 2 {uz+\)^ wnere Z2 - 1S de- 
termined from the scaling of the minimal gap along the 
path with respect to k [cf. Eq. (4) in Ref. 23, with a = 1 
and g?2 = 0]. Thus, Z2 = z in the quench across QCP A, 
leading to KZS, whereas Zi = 3 ^ z in the quench across 
MCP B, leading to anomalous scaling n ex (tf) ~ r" 1 / 6 . 
Given the above thermal initial condition, starting from 
Eq. (15) for the relative excitation probability, one finds: 

AP fe (t)=tanh(^^)|a , fe W| 2 , (29) 

where for both paths we simply write £k{h c ,Jc) to mean 
that critical parameter values are assumed at t = to. 
When T < e fe (/i c , 7c ), tanh ( tk{h £? a) ) w 1 and AP k (t) 
is the same as starting from the ground state of mode 
k. Thus, in order for the same ground-state scaling (ei- 
ther KZS or r -1 / 6 ) to emerge in the low-temperature 
limit, the condition T < efc(/i c ,7 c ) needs to be satis- 
fied for all the relevant modes. Since efc e (/i C )7c) = 0, 
this means that in the thermodynamic limit, the only 
allowed initial temperature is T = if a thermal state 
of H(h c ,j c ) is considered. In the opposite limit of high 
temperature, where T 3> efc(/i c ,7c), tanh ( ^^P^ ) 
ek{h c , 7 C )/(2T) ~ (k — k c ) z /T for modes k near k c . Upon 
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FIG. 9: (Color online) Exact scaling behavior of An ex (£/) in 
a quadratic adiabatic quench h(t) = 1 — j(t) = 1 — (t/r) 2 , 
starting from a thermal state at MCP B (to = t c = 0) toward 
the FM phase. The initial temperature T — 1000, yielding 
a linear fitting slope —0.663 ± 0.002, in good agreement with 
the value 2/3 predicted by Eq. (31). For comparison, the case 
of a ground-state quench is reproduced in the inset, with a 
linear fitting slope of —0.2190 ± 0.0006, which is also in good 
agreement with the predicted 2/9 exponent 23 . The data for 
different sizes (N = 800 and N = 1600) coincide up to 10~ 13 . 



integrating over the relevant modes and recalling Eq. 
(19), the relative excitation density is then: 



An ex (i/) 



1 



T 



AP k (t f )d d k 



-(d+z)vz/[z 2 (vz+l)} 



k z d d k 
(30) 



For the standard QCP A, this yields An ex (£/) ~ t~ 1 /T, 
recovering the result of Eq. (27), whilst An cx (t/) ~ 
t~ 1 / 2 /T in the multicritical quench across QCP B. In 
Ref. 23, we argued that the time-dependent excitation 
process in ground-state quenches need not be dominated 
by the critical mode k c for certain paths across MCPs 
and Pk = APk ~ k d ' 2 , with c?2 playing the role of an "ef- 
fective dimensionality exponent". For a thermal quench, 
it is interesting to note that, formally, one may inter- 
pret g?2 = z ^ in the above equation, also implying 
that the dominant contribution does not originate from 
modes around k c . In the high temperature limit, Pk{to) 
is, indeed, almost fully mixed for modes near k c , caus- 
ing the contribution of poo,k and pn.k to Eq. (15) to be 
nearly cancelled, and consistently leading to APfc(i) ps 
for those modes. 

The scaling prediction in Eq. (30) can be further gen- 
eralized to a non-linear thermal quench, whereby for in- 
stance h (t) = 1 — 7(t) = 1 — {t/r) a in the case of a quench 
away from the MCP B. When T = 0, Eq. (4) in Ref. 23 
yields 59 n ex (i/) ~ T -davz/[x 3 {auz+i)] _ Correspondingly, in 



FIG. 10: (Color online) Main panel: Scaling exponent 
AXX(tf) as a function of temperature T in a linear quench 
h(t) = 1 - j(t) = 1 - t/r away from the MCP B, starting 
with a thermal equilibrium state of H(h c ,y c ). Inset: Scaling 
exponent of AXX(tf) as a function of temperature in a linear 
quench h(t) — 1 + t/r away from the regular QCP A, starting 
with a thermal equilibrium state of H(h c = 1). In both cases, 
the system size TV = 800. 



the high-temperature limit, 



An cx (t f ) ~ —t 



-{d-\-z)a.vz f[z2 (avz-\-l)] 



(31) 



Exact numerical results for a quadratic quench (a = 2) 
are reported in Fig. 9, the inset corresponding to the 
ground-state T = case. Within numerical accuracy, 
the observed behavior is in excellent agreement with the 
predicted scaling, t~ 2 / 9 for T = and r~ 2 / 3 for high-T, 
respectively. 

We further examine how dynamical scaling is detected 
by other observables and how it is influenced by temper- 
ature away from the limiting regimes discussed above by 
considering the behavior of the spin correlator, AXX(t), 
defined in Eq. (18). Since XX does not commute with 
the Hamiltonian in Eq. (2), no analytical treatment is 
possible. Exact numerical results are presented in Fig. 10 
for both the regular and the multicritical QCPs A and B 
(inset vs. main panel, respectively), starting from the 
same thermal initial condition at criticality as consid- 
ered above. As the data show, similar features emerge 
in both cases: the scaling exponent of AXX(t), which 
is expected to be the same as for An ex , deviates from 
its zero-temperature value (—1/2 or —1/6, respectively) 
as soon as the temperature is nonzero, and as the lat- 
ter is gradually increased, it continuously changes until 
for sufficiently high temperature (T 3> €k(h c ,j c ), for all 
k € Kfj), it stabilizes at the value predicted by Eq. (31) 
(—1 or —1/2, respectively). All these observations are 
consistent with the predictions in the previous paragraph. 
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In summary, ground-state dynamical scaling (and KZS 
in particular) is fragile with respect to temperature fluc- 
tuations if the initial state is a thermal equilibrium state 
at criticality. In this case, the two situations where 
scaling exists are the zero-temperature and the high- 
temperature limit, with Eq. (31) holding in the latter 
regime. This requires all the relevant modes to either 
stay in their ground state or be highly mixed at the ini- 
tial time, which is a stronger condition in comparison to 
the ones identified in the previous sections for coherently- 
prepared (pure) excited states. From a practical stand- 
point, the high-temperature regime could potentially be 
relevant to liquid-state NMR simulators 50 . In order for 
tests of dynamical scaling/KZS in the low-temperature 
regime to be experimentally viable, however, the initial 
thermal state needs to be (or be prepared) sufficiently 
far away from criticality (e.g., \ho — h c \ 3> 1 for QCP 
A), in such a way that the condition T < A(k,ho) for all 
k £ Kpi can still be fulfilled with a non-zero temperature. 



B. Sudden quench dynamics and thermalization 

Sudden quenches have recently attracted considerable 
interest as a setting for probing the long-time dynam- 
ics of isolated many-body systems and the approach to 
equilibrium 36,38-42 . Since the quadratic Hamiltonian in 
Eq. (2) describes a simple (non-interacting) integrable 
model, it is well known that no thermalization can oc- 
cur in a proper sense, that is, the behavior of generic 
observables is not governed by a conventional statistical 
equilibrium ensemble 37,43 . The above investigations have 
nevertheless shown that information about the asymp- 
totic behavior of an appropriate subset of observables 
may still be encoded in a finite effective temperature T e g , 
independent on the fine details of the initial state and the 
dynamics but only determined by the total energy of the 
process. Let p(to) = Po an d Hj denote, respectively, the 
density operator describing the initial state of the sys- 
tem, and the final Hamiltonian after the (instantaneous) 
quench. Following Rossini et a/. 38 , the effective tempera- 
ture is defined by the requirement that the average energy 
of the initial state relative to the quenched Hamiltonian 
equals the one corresponding to a fictitious thermal state 
at temperature T c ff in the canonical ensemble, that is, 



Tr[p #/] = Tr[p Tog H f }. 



(32) 



Under the assumption that T = initially [that is, 
ground-state initialization in Eq. (32)], the emergence of 
effective thermal behavior has been related to the locality 
properties of different physical observables relative to the 
quasi-particle language that diagonalizes the model 39,40 . 
For a generic quench in a Ising chain, only non-local ob- 
servables (such as the two-point correlation functions of 
the order parameter) have been found to thermalize, with 
both their asymptotic average value and the finite-time 
transient being determined by equilibrium statistical me- 
chanics at T e ff. Remarkably, however, thermal behav- 
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FIG. 11: (Color online) Comparison between the long-time 
average quasiparticle excitation following a sudden quench 
ho M> hf starting from a thermal initial state at temperature 
T (dashed red) and the equilibrium value predicted by a ficti- 
tious thermal canonical ensemble at T e g (solid blue). Panels 
(a), (c), (e): sudden quenches to hf — h c — 1, hf — 1.01, 
hf — 1.5, respectively, with initial temperature T — 1.0. The 
behavior for a ground-state quench (T = 0, data not shown) 
is qualitatively similar, with deviations from the thermal pre- 
diction being further pronounced. Panels (b), (d), (f): sudden 
quenches to hf — h c = 1, hf = 1.01, hf = 1.5, respectively, 
with initial temperature T = 10.0. In all cases, N = 800, and 
the value of T e ff obtained from is Eq. (32) is also given. 



ior has also been established for certain local observables 
(the transverse magnetization per site, l/N^ja^, an< ^ 
the kink density, J\f) in quenches towards criticality, the 
long-time value being still univocally determined by T e g. 

Physically, it is clear that the concept of an effective 
temperature has a restricted validity and, for the model 
under investigation, it does not imply that an actual ther- 
mal ensemble emerges as a result of a sudden quench fol- 
lowed by free evolution under the quenched Hamiltonian. 
With that in mind, we further explore in what follows 
the emergence of effective thermal behavior in critical 
quenches, by focusing on a different local observable and 
by extending the analysis in two directions: first, initial- 
ization in a thermal state at finite T > and, second, 
sudden quenches to a multicritical QCP. 

Let us first consider a sudden quench of the magnetic 
field ho i— > hf in the Ising chain (7 = 1), starting from 
an initial state of the form given in Eq. (28), and focus 
on the long-time behavior of the number of quasiparticle 
excitations with momentum k. Since the corresponding 
observable commutes with the time-dependent Hamilto- 
nian, the long-time expectation value coincides 
with the one right after the quench. In order for the lat- 
ter to be consistent with the equilibrium value at T e ff, 
the following identity must hold: 

G°oo,fc(*o) - Pn,fe(io))|ao,fe| 2 + pu,k(to) + P33,k(to) 

= (1 + e +«*(*/.*r= 1 )/ 3 W)- 1 ) (33) 

where | o-o fe | 2 is the excitation probability of mode k due 
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FIG. 12: (Color online) Difference between the long-time 
quasiparticle excitation of the critical mode k c from its ther- 
mal equilibrium prediction as a function of system size for a 
sudden magnetic-field quench to h c in the Ising chain. An 
initial thermal state with temperature T = 1.0 is considered. 
The linear fitting slope is -0.99992 ± 3 x 10~ 5 . 
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FIG. 13: (Color online) Comparison between the long-time 
quasiparticle excitation following a sudden quench ho n* hf — 
h c = 1, 70 | — ► 7/ = 7c = towards the MCP B (blue) and the 
equilibrium value predicted by a fictitious thermal ensemble 
at T e g (red). The system is initially in a thermal state with 
temperature T — 10. Left: Initial state is the thermal state 
at ho = 2.0,7o = 1.0 (inside the PM phase). Right: Ini- 
tial state is the thermal state at ho = 0.0, 70 = —1.0 (inside 
the FM phase). Notice that due to the fact that the excita- 
tion probability of low-energy modes exceeds 1/2, T e g is much 
higher than in any other situation with the same initial T, cf. 
Fig. ll(b, d, f) and Fig. 13(a). 



to the quench and Eq. (11) has been used in the left hand- 
side. The right hand-side is the fermionic thermal equi- 
librium prediction Tr[py 7fc7fe]- Exact numerical results 
are presented in Fig. 11. Altogether, these data indicate 
that similar to the behavior of other local observables in 
a ground-state quench 39,40 , no effective thermalization is 
observed outside criticality [panels (c)-(f)], as expected. 
Even for a quench toward QCP A, however, the initial 
temperature T must be sufficiently high in order for our 
chosen observable to thermalize [panel (a) vs. (b)]. 

In order to gain physical insight into what distinguishes 
a critical vs. non-critical quench in our case, and un- 
derstand why effective thermal behavior fails to emerge 
outside criticality even for high initial temperature, it is 
useful to take a closer look at Fig. 11(d): clearly, the 
main difference between the equilibrium and the actual 
quasiparticle distribution arises from momentum modes 
close to k c . On the one hand, since A(fc c , hf) is the small- 
est gap at hf, the maximum quasiparticle excitation is 
expected to occur at k c from the equilibrium prediction 
[right hand side of Eq. (33)]. On the other hand, the peak 
of the observed distribution is located at modes close to 
k c , but not exactly at k c . Because the system is far from 
h c , note that the difference of poo.fc, Pn,k, and P33 t k f°r 
modes close to k c is negligible. Thus, the main difference 
is due to |ao,fc| 2 , which, as remarked, is the excitation 
probability of mode k at T = after a sudden quench to 
h c . Upon re-interpreting |ao.fc| 2 f»l - co.fc| 2 j it is possi- 
ble to make contact with the results shown in Fig. 8: 
clearly, the excitation probability of mode k c changes 
dramatically for ho close to h c , which suggests that k c 
does not contribute appreciably unless hf = h c . Instead, 
other modes close to k c can be excited for values hf w h c 
at which k c is not yet excited. Since the excitation con- 
tribution from such "quasi-critical modes" would then be 



larger than the one from k c . Eq. (33) would not hold. Ac- 
cordingly, the only way to enforce the validity of Eq. (33) 
is through a sudden quench towards h Cl as observed. 

Having clarified why criticality is essential, we turn to 
assess whether the requirement of a sufficiently high ini- 
tial T may be related to the finite system size or will 
persist in the thermodynamic limit. We focus on a sud- 
den quench ho = 3.0 1— > h c at T = 1.0, and analyze 
how the long-time average of the total quasiparticle den- 
sity l/Nj^k 7fc7fc deviates from the thermal equilibrium 
prediction at T e g as N is increased. While we find that 
the observed deviations are practically constant over the 
range of N explored (data not shown) , the difference be- 
tween (7I 7fc c ) and its corresponding thermal prediction 
at T e ff does decrease with increasing N: as seen in Fig. 12, 
such a difference Ajl o 7k c ~ AT" " 992 at T = 1.0, im- 
plying a vanishing difference and effective thermal be- 
havior also at low temperature for the critical mode as 
A" — > 00. This property, however, stems from the fact 
that the gap of k c closes in the thermodynamic limit, 
which is not true for the gap of other modes. For either 
the number of quasiparticles in a generic mode or for the 
total quasiparticle density, we thus conjecture that even 
in the thermodynamic limit, thermal behavior will be ob- 
served following sudden quenches to the QCP A only if 
T ^> Afe(ft, c ,7 = 1) for all the relevant modes. 

In view of the peculiar features that distinguish a mul- 
ticritical QCP, as reflected in particular in anomalous 
scaling behavior 23 , it is not obvious whether the above 
condition would still suffice for the same observables to 
thermalize in a sudden quench towards MCP B. Exact 
results for two sudden multicritical quenches of the form 
{h Q = 1 + 70 (->• hf = 1 + 7/,7o 7/) are given in 
Fig. 13, starting from a thermal state at high tempera- 
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ture: specifically, MCP B is both reached via a sudden 
quench from the PM phase (left panel) and via a sud- 
den quench from the FM phase (right panel). Contrary 
to the high-temperature scenario for the regular QCP A 
[Fig. 11(b)], no thermal behavior emerges, the observed 
expectation value (T k Jk) for modes close to k c being sig- 
nificantly smaller or larger than the thermal equilibrium 
prediction, respectively. 

This anomalous long-time behavior can be traced back 
to the asymmetry of the impulse region along the con- 
trol path, as sketched in Fig. 1 (bottom). Following 
Ref. 23, the location of the minimum gap for each mode 
k along the path h = 1 + 7 is determined by requiring 
9A fc ( 7 ,l + 7)/9 7 = 0, that is, 

y(k) = (cos k - 1) /(l + sin k 2 ) < 0, 

which indicates that the center of the impulse region is 
largely shifted into the FM phase for each k. As a result, 
after a sudden quench to the MCP B from the FM phase, 
the excitation probability of low-energy modes tends to 
be enhanced above 1/2, whereas for a sudden quench to 
MCP B from the PM phase, the excitation probability 
of low-energy modes tends to be suppressed below 1/2. 
Since the thermal equilibrium value is close to 1/2 in 
the high-temperature limit for low-energy modes, ther- 
mal behavior is not realized in either quench process. 

Based on the above results, we conjecture that quench- 
ing toward the center of the impulse region is a necessary 
requirement for 7^7^ or l/-^5D fc 7^7fc to thermalize fol- 
lowing a sudden quench. While typically this is the case 
in a quench to a regular QCP (for instance, a sudden 
quench of h to QCP A at fixed 7 = 1), for a sudden 
quench to MCP B along the path h = 1 + 7, the loca- 
tion of the minimum gap (hence the center of the impulse 
region) is different for each mode k, preventing thermal- 
ization to be possible along this path irrespective of the 
final values hf , 7/ . More generally, we expect the above 
requirement to be necessary for local observables other 
than those examined here. While this goes beyond our 
current scope, it would be interesting to verify, for in- 
stance, whether the transverse magnetization or the den- 
sity of kinks would still effectively thermalize in a multi- 
critical quench to QCP B from the ground state. 

We also remark that in a recent work 41 , general conclu- 
sions have been reached for the equilibrium distribution 
after a sudden quench, predicting, in particular, effec- 
tive thermal behavior for generic observables when the 
quench is performed around a non-critical point, and 
poor equilibration otherwise. While at first these re- 
sults seem to contradict both our present conclusions 
for critical quenches towards QCP A in the appropri- 
ate temperature regime as well as earlier results for zero 
temperature 39,40 , a crucial assumption in Ref. 41 is a 
small quench amplitude, causing only a small number of 
excited states to effectively contribute around a QCP. 
The opposite condition is implied throughout our dis- 
cussion, the sudden quench amplitude being in fact large 



enough for the number of excited states involved in a crit- 
ical quench to outweigh those involved in a non-critical 
one (cf. Fig. 8). In the light of that, we also conjecture 
that having a sufficiently large number of states involved 
in the excitation process is a general necessary condition 
for effective thermalization after a sudden quench. 



V. CONCLUSIONS 

In summary, we have addressed how different aspects 
of many-body non-equilibrium dynamics depend upon 
initialization in a state other than the ground state for 
a class of one-dimensional exactly solvable XY models. 
Our main findings may be itemized as follows: 

• Dynamical scaling: initial pure excited states. 
Provided that the non-equilibrium response of the sys- 
tem is characterized in terms of suitable relative indica- 
tors (such as the relative excitation density), adiabatic 
quench dynamics can still encode the ground-state equi- 
librium critical exponents for a large class of initial energy 
eigenstates as well as for pure excited states prepared by 
a sudden parameter quench. A crucial role is played by 
how the initial excitation is distributed over the set of 
relevant quasi-particle modes that effectively evolve in an 
adiabatic quench. In particular, a unifying criterion that 
ensures the emergence of KZS in both the above scenar- 
ios in the thermodynamic limit is obtained by requiring 
that the majority of the relevant modes share a common 
initial excitation pattern, as expressed by Eq. (22). 

Our results recover ground-state scaling when no exci- 
tation is initially present, but they also allow for the crit- 
ical exponents of the ground-state QPT to be encoded in 
the scaling behavior for highly energetic initial configu- 
rations, where most of the relevant modes are fully or 
half-excited. While this makes contact with similar con- 
clusions on critical entanglement scaling in excited en- 
ergy eigenstates recently obtained for static QPTs 53 , it 
confirms that only the set of relevant modes is important 
in dynamical scenarios, as opposed to the full static mode 
set. Beside being supported by exact numerical methods 
and analytical derivations in limiting regimes, a justifica- 
tion of the proposed criterion has also been obtained for 
the case of excited-eigenstate initialization by suitably 
extending the perturbative (first-order) AR approach we 
previously employed for ground-state continuous QPTs. 

• Dynamical scaling: initial mixed states. In 
general, more restrictive conditions on the distribution of 
the initial excitation over relevant modes must be obeyed 
for universal dynamical scaling to emerge in adiabatic 
quench dynamics that originates from a statistical (inco- 
herent) mixture as compared to a (coherently prepared) 
pure state. In particular, two distinct scaling regimes 
have been identified for an initial thermal ensemble at 
a finite temperature T, depending on whether the latter 
is very low or very high relative to the relevant quasi- 
particle energy scale, and leading to KZS t -1 / 2 vs. r _1 
for a standard QCP, respectively. Since in both cases 
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all the relevant modes must share a common excitation 
pattern if the initial thermal state is prepared at crit- 
icality, this implies that KZS is fragile against thermal 
fluctuations in this case, the scaling exponent deviating 
from the KZ prediction as soon as T ^ 0. From a prac- 
tical standpoint, it is however important to note that a 
finite range of temperatures can still support KZS if the 
system is initially at thermal equilibrium sufficiently far 
from criticality. General predictions for scaling behavior 
in adiabatic thermal quenches involving a MCP have also 
been obtained [cf. Eq. (31)], and verified to be consistent 
with exact numerical results. 

• Effective thermalization. Effective thermal be- 
havior may emerge in the relaxation dynamics of the 
quasiparticle density following a sudden quench from a 
thermal state under appropriate conditions. Specifically, 
the long-time expectation value of this observable is de- 
termined by a fictitious thermal equilibrium ensemble at 
temperature T ft provided that i) the system is quenched 
toward the center of the impulse region, and ii) the ini- 
tial temperature is sufficiently high with respect to all the 
relevant gaps. For a standard QCP, the first requirement 
is met by a sudden quench toward criticality, which has 
been found sufficient for local observables such as the 
transverse magnetization per site and the kink density 
to thermalize starting from the ground state 39,40 . Our 
results indicate that, in general, condition i) alone need 
not suffice for arbitrary local observables. While require- 
ment ii) may be taken to be in line with what expected for 
a free (integrable) theory 36,43 , it remains an interesting 
open question to precisely characterize what subclass of 
local observables may ehxibit effective thermal behavior 



under the sole condition i) 60 . 

Our results additionally show that for certain ob- 
servables (such as the quasiparticle density), effective 
thermalization may fail to occur altogether (or possi- 
bly require yet more stringent requirements) for sud- 
den quenches to a multicritical QCP. Physically, we have 
traced this behavior back to the existence of quasicriti- 
cal (path-dependent) energy states and the correspond- 
ing shift of the impulse region, which also underlies the 
emergence of anomalous scaling exponents 23 . In this con- 
text, an interesting next step would be to examine the 
thermalization behavior of other local observables as con- 
sidered in Ref. 39,40. 

While the above analysis provides a more complete 
picture of non-equilibrium dynamics in a paradigmatic 
class of spin chains than available thus far, it remains 
a main open question to understand how crucially our 
results rely on the XY chain being an exactly solvable 
non-interacting model. From this point of view, it would 
be worthwhile to explore, for instance, whether dynami- 
cal critical scaling may still exist for finite-energy initial 
states in non-integrable models, or even in more com- 
plex but still integrable systems such as a Bethe-Ansatz 
solvable one-dimensional Heisenberg XXZ chain 53 or an 
infinitely coordinated Lipkin-Meshkov-Glick model 30 . 
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